Fourth moments reveal the negativity of the Wigner function 
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The presence of unique quantum correlations is the core of quantum information processing and 
general quantum theory. We address the fundamental question of how quantum correlations of 
a generic quantum system can be probed using correlation functions defined for quasiprobability 
distributions. In particular we discuss the possibility of probing the negativity of a quasiprobability 
by comparing moments of the Wigner function. We show that one must take at least the fourth 
moments to find the negativity in general and the eighth moments for states with a rotationally 
invariant Wigner function. 



I. INTRODUCTION 

The von Neumann projective detection scheme [l[ al- 
lows to measure simultaneously only commuting quan- 
tum observables. It is impossible to define a useful prob- 
ability of measurement outcomes for noncommuting ob- 
servables, like position and momentum, consistent with a 
projective measurement. To circumvent this problem, it 
is possible to define a quasiprobability - the Wigner func- 
tion (2j, consistent with the projective measurement but 
taking sometimes negative values. Because of the Heisen- 
berg uncertainty principle, the Wigner function cannot 
be directly measurable in the von Neumann sense. How- 
ever, one can find it by indirect measurements like the 
homodyne tomography [3j, assuming a special detection 
scheme and using the inverse Radon or Abel transfor- 
mations. In this scheme it is essential to use efficient 
photon counters. The negativity of the Wigner function, 
with the position and momentum replaced by the electric 
and magnetic field, has been demonstrated using single 
photons and other states 0, [5[ . 

Recently, it has become possible to measure the mo- 
ments of the microwave electromagnetic field @, @] ■ The 
microwave frequencies make it difficult to build efficient 
photon counters although there exists a proposal [8] . Us- 
ing instead linear amplifiers leads to a large detection 
noise and, contrary to the homodyne detection, one can- 
not reconstruct the complete Wigner function. Hence, it 
is relevant to ask if we can show the negativity of the 
Wigner function having access only to correlations of po- 
sition and momentum. 

The detection schemes based on linear amplifiers are 
easier to interpret in terms of the quantum weak mea- 
surement 0] than in the traditional projective scheme. 
The interpretation of weak measurements leads to para- 
doxes similar to the negative Wigner function. In the 
weak measurement picture, the detector just adds a large 
white noise of the outcome. The detection noise contains 
no information about the measured system and can be 
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subtracted by a deconvolution of the resulting probabil- 
ity distribution. However, the remaining distribution is 
a quasiprobability, taking sometimes negative values. In 
the case of momentum and position the quasiprobability 
corresponds to the Wigner function [RM12l |. The decon- 
volution is in principle experimentally feasible but since 
usually one has access only to correlation functions, it 
is useful to develop schemes for them to detect quan- 
tum correlations. Generally we have shown previously 
that the knowledge of a finite number of moments of a 
quasiprobability can show its negativity [l2j . 



In this paper, we discuss how the negativity of the 
Wigner function can appear as violation of an inequal- 
ity involving only finite moments. We emphasize that 
there exist examples violating the Cauchy-Bunyakovsky- 
Schwarz inequality by fourth or higher moments for 
position and momentum (or second moments for their 
squares or intensities) 1 131.1 141 but only using the Glauber- 
Sudarshan P-function [15| , whose negativity is much eas- 
ier to demonstrate, compared to Wigner function. We 
will show in the following that it is possible to prove the 
negativity of the Wigner function using only the fourth 
moments for position and momentum although the vi- 
olating state is complicated. For rotationally invariant 
states, mixtures of Fock states for a harmonic Hamilto- 
nian with a rotationally invariant Wigner function, we 
show that one needs at least the 8 th moment to show the 
negativity of quasiprobability. Additionally, we provide 
a connection between continuous time-resolved measure- 
ment of position and momentum which permits to find 
moments of a single-time Wigner function by deconvolu- 
tion in a harmonic system. 



The paper is organized as follows. We start from re- 
calling basic properties of Wigner function, including its 
negativity. Next, we demonstrate violation of classical in- 
equalities by fourth moment for a special state and eighth 
moments for a rotationally invariant state. Finally, we 
present a way to find moments in continuous measure- 
ment by deconvolution and close with conclusions. 
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II. WIGNER FUNCTION 

We begin with the standard construction of the Wigner 
function. Let us define the dimensionless Hermitian po- 
sition and momentum operators, x, p with the property 
xp — px = [x,p] = i (h is incorporated in the definition 
of x and p). One can always make a rescaling x — > X x x, 
p — > X p p later to recover the proper dimension. Instead 
of x and p, in optical systems one can use the electric and 
magnetic fields, E and B for a given mode. The Wigner 
function is defined as 0, 

w(x, P ) = f ^| e - i (x»+ep) Tre i(xi+ep) /5 . 

J (271") 

Here p is the density matrix of the system. One essen- 
tial property of the Wigner function is that its marginals 
correspond to the standard projective probability, 

p(x)=TrpS(x-x) = jw( X> p)d P) 

pip) = TrpS(p-p) = J W(x,p)dx. (2) 

Similar relations can be found for arbitrary linear com- 
binations of x and p, used e.g. to experimentally extract 
the Wigner function We stress that, contrary to the 
P— function, the Wigner function definition does not in- 
volve any characteristic harmonic frequency u). Neverthe- 
less it is convenient to represent p in terms of coherent 
or Fock states defined for a harmonic oscillator. 

The Wigner function can never be measured directly, 
which follows from the Heisenberg uncertainty principle 
and the consistency of quantum measurements. It cannot 
be measurable as a probability because it takes negative 
values. To show this we introduce recursively the Fock 
basis states \n + 1) = aJ\ri)/\/n + 1 for n — 0,1,2,... 
and (0|0) = 1, d|0) = 0. They form the complete set of 
eigenstates of the Hamiltonian 

H = hwa^d 7 a = (x + ip)j\pl, (3) 

with the property [a, fit] = 1. In this notation we have 
set to = lu^ 1 . We have H\n) = ftwn\n) and 

W{x,p) = e-( x2+p2 \2{x 2 + p 2 ) - l)/ir (4) 

with W(0,Q) = -1/tt < for p = |1)(1|. In general, 
for pure states p = the Wigner function is al- 

ways somewhere negative unless the state has a Gaus- 
sian wavefunction [161 ]. The directly measurable quan- 
tity is W(x(t)) which allows to retrieve the full Wigner 
function by a numerically difficult inverse Radon trans- 
form (or Abel transform if the state can be assumed to 
be rotationally invariant) Q- The statistics needs a lot 
of repetitions of the experiment. Moreover, it needs a 
strong, projective measurement, excluding any external 
noise. This condition is possible to fulfill theoretically but 
the experimental implementation relies much on a hardly 



verifiable detection scheme (based on efficient photon- 
counters). More common is the situation with some, 
in many cases even large detection noise, which on the 
other hand can be well estimated and subtracted Q as 
we show in Sec. V. This latter method makes it impossi- 
ble, however, to find the whole Wigner function. In this 
scheme one can only find a finite number of moments of 
the Wigner function. We stress that the moments of the 
Wigner function are naturally measured in a continuous 
weak measurement (and not the P- and the Q-function) . 
Furthermore, as we show in Sec. V by a suitable process- 
ing of the continuous signal of position only, all moments 
involving momentum can be extracted as well. 

Motivated by the continuous measuring scheme it is 
natural to pose the following questions. Can we say some- 
thing about the negativity of the Wigner function from 
the knowledge of a finite number of its moments? What 
is the smallest order of the moment to show that the 
Wigner function is negative? How complicated must the 
violating state be? Below we answer all these questions. 

III. NONCLASSICAL MOMENTS - GENERAL 
CASE 

An arbitrary probability distribution must satisfy the 
generic inequality 

(f 2 (x, P )) >0. (5) 

Once the above inequality is violated, the candidate for 
a probability must be somewhere negative and fails the 
test to be a real, positive probability. In the case of the 
Wigner function we consider the following expansion 

f(x,p) = ^ c nmX n p m = C + C W X + C m p + 

I vm 

C2 x 2 + cuxp + c 02 p 2 + c 3Q x 3 + ... (6) 

with real c nm for n, m > 0. We say that / is of the order 
k if there exits a nonzero c nm for 77 + 777 = k and all 

Cnm — for 77 + 777 > k. 

We focus on the following task: find the state p that 
violates 

(f 2 (x,P))w>0 (7) 

for the / of the lowest possible order. One can easily 
show that k > 1. If f(x,p) = cq + ciqx + coip then 

(f{x,p))w =Tr/ 2 (i,p)p>0. (8) 

Next, the case k = 2 can be reduced by linear transfor- 
mations on the phase space, namely rotating (x,p) — > 
(xcosip + p sin 0, p cos — xsmifi), squeezing (x,p) — > 
(gx,g~ 1 p) and shifting (x,p) — > (x + xq,p + pq), which is 
justified in both the classical and the quantum case, to 
f a (x,p) — x 2 +p 2 + c or fb{x,p) = 2xp + c . Other cases 
either reduce to f a ,b or will never violate the inequality 



3 



(fTj) . It is enough to consider pure states p = \tp)(tp\ be- 
cause violating a mixed state must contain a violating 
pure state in its decomposition. Remember that choos- 
ing the Fock space of the harmonic oscillator we do not 
lose generality, since an arbitrary state can be written in 
terms of Fock states even if the Hamiltonian is different 
from 



For f a , one can show that x 



2a*a and 



(2a*a) w = Tr(2n + l)p, 

((2a* a) 2 ) w = Tr(4h 2 + 4n + 2)p, (9) 

with h — a^d. Hence, we can consider only the |n) states 
but An 2 + An + 2 > (2n + l) 2 so /„ will never violate (J7J. 
We have 2xp — i(a* 2 — a 2 ) and 



((a* 2 -a 2 )) w = Tr(a+ 2 - a 2 )p, 



(10) 



((o* 2 - o 2 ) 2 )w = Tr [(a 4 + a t4 )p] -2((a*a) 2 ) w . 
This gives also the bound (f 2 )w _ —1 following from 

(fb)w + 1 = Trp(i(aV - a 2 ) + c ) 2 > 0. (11) 

We can consider separately l^) made of odd states |2fc+l) 
and of even states \2k) because fb preserves parity, since 
operators appear pairwise. 

The search for the violation of ([7]) reduces to the di- 
agonalization of a quadratic form, parametrized by Co, 
which is easily accomplished numerically. The outcome 
of the search is that we need at least 5 states of Fock 
space and the largest violation of ([7]) occurs for Co = 0. 
The simplest state violating (J7J with /(, has the structure 

\tl>) = «o|0) + v A \A) + v 8 \8) + v 12 \l2) + v w \m). (12) 

Since Co = and xp shifts the state number by 2, we 
omit states |4fc + 2). The coefficients (vq, Vi, v$, vx 2 , v\%) 
form the eigenvector of the matrix 



/4 ft \ 

B x A x B 2 ' 

B 2 A 2 B 3 

B 3 A 3 B 4 

V B 4 A 4 / 



with the smallest eigenvalue for 

A k = 32fc 2 + 8fc + l, 

B k = -y/4k(Ak - l)(4fc - 2)(Ak - 3). 



(13) 



(14) 



The fact that the matrix has a negative eigenvalue follows 
from the negativity of its determinant, equal —10447775. 
The numerical diagonalization gives X m in = ((^xp) 2 )w = 
—0.036 and its normalized eigenvector is given by 
(v ,v A7 v 8 ,vx 2 ,vi e ) = (0.973,0.206,0.0897,0.042,0.0161). 
Going to more complicated states allows for an even 
stronger violation. However, we stress that our numerical 
search has shown, that at least a superposition of 5 Fock 
states is necessary. Hence, we draw the important con- 
clusion, that showing the non-classicality of the Wigner 



function using only moments requires at least fourth mo- 
ments and states involving coherent superposition of at 
least 5 states. 

The violation for a general state is bounded by 
Eq. (fTT|) . It is interesting to note, that one can find a 
state with (f 2 )w arbitrarily close to —1. We can con- 
veniently describe such a state by its wavefunction in 
position space \ip) = J dx^(x)\x), with (x\y) = 5{x — y), 
x = x and p = —id/dx. We write (llip in the form 



(fb)w + l= / \Hx)\ 2 dx 



(15) 



and 



4>(x) = (xp+px+co)ijj(x) = (co — i—2ixd/dx)i()(x). (16) 

To get the minimum of (f 2 )w we should make 4>(x) as 
small as possible. However, trying to get (f 2 )w — — 1 
implies = and ip(x) ~ |a;| _ ( 1 + lCo )/ 2 , which cannot be 
normalized. Instead, we take a regularized wave function 



^(x) = e-l*l/ 2 |a:| e " 1/2 /v^r(2i) 



(17) 



with e > and the Gamma function T(z). The Wigner 
function reads 



W(x,p) 



2 c 



Ig/gj 
Trr(e) 



nJ e (2xp) 2K t (\x\ +2ixp) 



Ixp^e 



!■'•! 



Re 



(ipx+ \x\/2Y 

(18) 

with W - Jo(2xp) - Y (2xp) for e - and J,K,Y de- 
noting Bessel functions. We get <f>(x) = —i(2e — \x\)ip(x), 
which gives ((2xp) 2 )w = 2e — 1. This result shows that 
we can get arbitrarily close to the limiting value — 1 with 
e — > 0, however, the case e = is impossible as the wave 
function would not be normalizable. We stress that pTj) 
represents only an example of a state with (f 2 )w — > — 1 
and one can certainly find also such states in the family 
of analytic functions, e.g. a sequence of states in finite 
Fock spaces. 



IV. NONCLASSICAL MOMENTS AT 
ROTATIONAL INVARIANCE 

We ask the further question: What is the lowest or- 
der of / that violates (UJ for rotationally invariant states 
(considering the Wigner function in phase space)? The 
rotational invariance is intimately connected with the 
harmonic dynamics ([3|), which induces a rotation fre- 
quency uj. It is important since it is often difficult to 
control experimentally the phase of the photonic state, 
corresponding to the polar angle in the phase space. If 
the phase is smeared out then the state is rotationally in- 
variant, described by a mixture of Fock states |n), namely 
p = ^2 n Pn\n)(n\. Again, it is enough to focus on pure 
Fock states. Next, the parts of / with different parity 
with respect to x or p can be considered separately as 
their cross products integrate out to zero for Fock states. 
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The function /& reduces to f a for Fock states so it cannot 
violate ([7]). The only nontrivial case with k = 3 is then 



f c {x,p) = x{x 2 + p 2 ) + (c 30 - 1).t 3 + c 10 s 



(19) 



One can find that for a rotationally invariant distribution 

(/c> = (5cio + 2 C 30 + l)(r 6 )/16 

+(3c 30 + l)ci (r 4 )/4 + c 2 (r 2 )/2, (20) 

where r 2 — x 2 +p 2 . The straightforward analysis of the 
above quadratic form gives its minimum located at 



C30 



cio 



l )(r 



ill) = 



3(r 4 ) 2 -2(r 6 )(r 2 ) 

10(r 6 )(r 2 ) -9(r 4 ) 2 ' " iu 10(r 6 ) (r 2 ) - 9(r 4 \ 2 
(r 6 )((r 6 )(r 2 )-(r 4 ) 2 ) 



2(10(r 6 )(r 2 ) -9(r 



(21) 



1 w 



On the hand, for the Fock state \n) one finds (r 2 
2n+l, (r 4 ) w = 4n 2 + 4n + 2 and (r 6 ) w = 8n 3 + 12n 2 + 
16n + 6 so (r 6 ) w (r 2 ) w - (r 4 ) 2 ^ = 12n 2 + 12n + 2 > 0. 
Hence, (/ 2 ) will be never negative for Fock states. 

The negativity appears for k = 4, namely for fd(x, p) = 
r 4 + c 2 or 2 + Co, taking the state |1). One gets 

(fd)w = 10c^ + c 2 + 6c 20 c + 84c 20 + 20c + 216. (22) 

The minimum of the above function is located at C20 = 
— 12, Co = 26 and equals (f3)w = —28. This is our sec- 
ond result: we need the eighth moments and not less to 
show the nonclassicality of rotationally invariant states. 



V. MOMENTS FROM CONTINUOUS 
MEASUREMENT 

The concept of the Wigner function is very general as 
it only relies of the fact that two observable do not com- 
mute. However, in many real experiment the observables 
are measured continuously and, hence, it necessary to 
extend the concept of a quasiprobability to continuous 
time-resolved measurements. To this end, we define a 
time- extended Wigner functional, in accord with general 
considerations [10l - fl2l ]. by a path integral 



W[x,p] = J D X D£ e 



S\x,S]-J idtlx(t)x(t)+£(t)p(t)] 



e S[x,£l = ^ r j- e S ^t[ x (t)x(t)+at)P(t)]/2 p x (23) 
je/«tt[x(t)*(*)+e(t)£(t)]/2 

where T(T) denotes (anti)time ordering, A(t) denotes 
an operator in the Heisenberg picture and p is the initial 
density matrix. Such a functional can be measured as a 
convolution with some detection noise [l2j 

p[x,p} = I Dx'Dp'p d [x-x',p-p']W[x',p'}, (24) 



where pd corresponds to the large detection white noise, 
independent of the state of the system. This corresponds 



to the experimental detection scheme, in which moments 
of W can be found by deconvolution 0, [llj ■ The mea- 
surable quantities are the cumulants 



({A 1 ---A n )) = (-i) 
and the moments 

(A 1 --.A n ) = (-i) n 



n d n S(xi,...,Xn) 



dxi--- dxn 

gn e S(xi,...,Xn) 



(25) 



(26) 



dxi--- dxn 

where the cumulant generating function (CGF) is defined 
S(xi,...,Xn)=ln<e i ^+- + ^). (27) 

Knowing the cumulants of p and pd one can find the 
cumulants and moments of W due to the relation between 
the CGFs 



S p — Sd + s 



w ■ 



(28) 



Now, let us focus on a harmonic oscillator described 
by the Hamiltonian ([3]). In the Heisenberg picture a(t) = 
e~ lut a. Due to the fact that all operators in (|23|) com- 
mute to c-numbers, one can now write 

e s[x,€\ — T]. e / idt[ x (t)x(t)+t(t)p(t)] p. (29) 

One can see also that the Wigner functional ([2^| vanishes 
outside classical trajectories, namely 

(da(t)/dt + iua(t))W[x,p] =0, a = (x + ip)/V2. (30) 

It leads to a simple picture: once the initial distribu- 
tion W(x,p) is known, the general functional VV[ar,p] is 
obtained by evolving along classical trajectories. There- 
fore, for a harmonic Hamiltonian, the complete Wigner 
functional is not necessary for the theoretical calcula- 
tions as it reduces to the single-time Wigner function. 
However, the classical evolution is helpful for the exper- 
imental detection schemes. It is worth to collect all data 
from a continuous (weak) monitoring of x(t) and/or p(t) 
and construct the initial Wigner function by a reverse 
mapping of the classical dynamics. 

Let us introduce the notation for the marginal distri- 
butions 



W[x] = J DpW[x,p], 



V\?{xi, . . . ,x n ,px, . . . ,p m ) = / DxDpW[x,p] x 



5(xi - x(ti)) ■ ■ ■ S(x n - x{t n )) x 
S(p 1 -p(t' 1 ))---S{ Pm -p(t' m )). 



(31) 



We shall switch the notation Xk — > x(tk) when it is un- 
ambiguous. It is easy to show that 

W[a:(0),a;(7r/2w)] = W[z(0),p(0)] = W(x,p), (32) 
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so one can obtain the usual Wigner function also from 
only position measurement, but at different times, as 
p(0) = x(tt/2uj) follows from the equation of motion (pO)) . 
More generally, W(x,p) can be obtained also from the 
appropriate Fourier components, namely 

W(x,p) = J Dx' sfx- 2cos(ujt)x'(t)/t^ x 

S (p-J 2sm(iot)x'(t)/t^W[x'}, (33) 

where to ^ is the total averaging time. We stress 
that the averaging time is limited by the backaction of 
the detector, which protects - by means of the detection 
noise - the Wigner function from its exact measurement 
in a single run of an experiment. 



but at least of fourth order. As also shown in other con- 
text we confirmed that fourth order correlations are nec- 
essary to show nonclassicality of quantum mechanics [HJ ■ 
We believe this is fundamentally important for probing 
nonclassicality by means of weak measurement, when one 
has access only to moments. We are convinced that our 
results will stimulate an effort to measure higher order 
correlations in optical cavities by means of linear detec- 
tors and amplifiers even if they add large detection noise. 
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